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Geometrical effects on spin injection: 3D spin drift diffusion model 
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We discuss a three-dimensional (3D) spin drift diffusion (SDD) model to inject spin from a ferromagnet (FM) 
to a normal metal (N) or semiconductor (SC). Using this model we investigate the problem of spin injection 
into isotropic materials like GaAs and study the effect of FM contact area and SC thickness on spin injection. 
We find that in order to achieve detectable spin injection a small contact area or thick SC samples are essential 
for direct contact spin injection devices. We investigate the use of thin metal films (Cu) proposed by S.B. 
Kumar et al. and show that they are an excellent substitute for tunnelling barriers (TB) in the regime of 
small contact area. Since most tunnelling barriers are prone to pinhole defects, we study the effect of pinholes 
in A10 tunnelling barriers and show that the reduction in the spin-injection ratio (7) is solely due to the 
effective area of the pinholes and there is no correlation between the number of pinholes and the spin injection 
ratio. 



I. INTRODUCTION 

Spintronics has been a rapidly growing field from the 
past few years. The main interest in spintronics has 
arisen because of a spin-interference device proposed by 
Datta and Das^, which is based on spin precession con- 
trolled by a gate voltage. Since then a lot of modified 
devices have been proposed, but in all these devices the 
main problem is of efficient spin injection from the ferro- 
magnet (FM) to the normal metal (N) or semiconductor 
(SC). In order to understand efficient spin injection it 
is necessary to understand the coupling between charge 
and spin currents, which was first described by Aronov 2 
and later developed by Johnson and Silsbee in terms of 
thermodynamic processes^. P.C. van Son et al4 later 
proposed a much simpler linear response model based 
on spin drift and diffusion (SDD) to describe transport 
across FM-N/SC interfaces. SDD model was successfully 
applied to current perpendicular to plane geometries of 
giant magneto-resistance by Valet and Fert and they also 
established the connection between the diffusive model 
and the Boltzmann equation^. The problem in injecting 
spin from FM to SC is due to the difference in the con- 
ductivities and spin diffusion lengths of the materials and 
is termed as the conductivity mismatch problem 6 , which 
was solved by Rashba using spin dependent boundary 
resistance^. 

Although SDD has been used to describe spin trans- 
port across many local^ and non-local geometries^—, 
it is applied as a ID theory. It was first shown by T. 
Kimura et al. that a ID theory is insufficient to describe 
a 3D experimental geometry^. A 2D extension of the 
SDD model by Ichimura et al. has been used to study 
the spatial variations of the spin current and the elec- 
trochemical potentiate and also a quasi 3D model based 
on spin dependent resistive elements (SDRE) was pro- 



posed by Hamrle et al., where each SDRE follows ID 
equations^. Although these models are extensions to 
the ID case they are still insufficient to take into account 
all the effects encountered in an experimental 3D device. 
In this paper we will first describe the SDD model in 
3D with the appropriate boundary conditions and point 
out the shortcomings of the ID model. We then apply 
our model to describe the effects of SC height and the 
FM contact area on the spin-injection efficiency of the 
device. Rere we will discuss key issues such as the possi- 
bility of direct contact spin injection device and the use 
of thin metal films as substitutes for tunnelling barriers. 
Our main findings are that when device parameters are 
smaller than the order of the spin diffusion length, the 
spin injection ratio is influenced dramatically which is 
not predicted by the simplistic ID model. Finally, we 
will discuss pinhole defects in tunnelling barriers and the 
effect of pinholes on the spin injection ratio of the device. 



II. THEORY 

In order to derive the spin drift diffusion equations we 
assume that far from the interface at temperatures lower 
than the Curie temperature most scattering events will 
conserve the spin direction and thus the spin up and spin 
down electrons will flow almost independently of each 
other—. Also if the spin scattering occurs at much longer 
time-scale than other electron scattering events we can 
define the electrochemical potentials fi t and /i^ for both 
the spin channels. Thus in the linear response regime 
the current carried by the spin-up (j t ) and spin-down 
(jj.) channel is given by Ohm's law: 



Jt,4- 



(1) 
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where, er t 4 = a(l ± a)/2 is the spin dependent electrical 
conductivity^ and e (> 0) is the electron charge. Near 
the interface the spin can diffuse from the up-spin channel 
to the down-spin channel and thus the coupling of the two 
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spin channels is given by the diffusion equation: 



(2) 



where, D is the Diffusion constant and V Dt = A, the 
spin diffusion length. 

Now in order to simplify notation we use the following 
transformations 8 



C = Mt - M-m 
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where T is the spin-injection ratio and n x is taken as the 
normal to the surface along the flow direction. 
Eq. JJJ) and © then transform into: 
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where, Act = er t — er 4 and J = j t + is the total current 
through the system. 

The boundary conditions for the transformed equa- 
tions are-^: 

C4,-C*|o = 2r„ fr-^J (J-nOe, (9) 



CU~ = 0, (10) 



AT, 

Z,\, - Z F \ = r„ ( 1 - — r ) (J • rW)e, (11) 
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(12) 



(13) 



E/(4E t E + ) is 

the effective contact resistance and n 2 is the normal to 
the boundary of the domain. Here subscript denotes 
the interface and f2 is the domain of the device^. In 
addition to Eq. ©, dTUJ) , ([IT]). (IT2l which are generally 
used for the ID model, we use the Eq. (fT3")) for the 3D 



FIG. 1. (Color Online) Spin injection ratio as a function of 
the bridging distance d for a three-terminal device. Square 
(red-online) represents the ID result. 



model. The additional boundary condition ensures that 
no spin current leaks out of the device. 

We solve Eq. © and using the program FF3D21, 
which employs fictitious domain finite element method 
(FEM). Fictitious domain method allows us to change 
the device geometry without actually changing the grid, 
hence in the case of defects in tunnelling barriers this 
turns out to be the most viable option. In order to ob- 
tain the spin injection ratio we solve the equations in 
the two regions (FM and SC/N) separately and itera- 
tively vary the boundary conditions till a convergence of 
1CP 6 in the spin injection ratio is achieved. To ensure 
that all the results are well converged with respect to the 
grid parameters we make sure that for boundaries where 
Eq. (fT0|) needs to be satisfied the boundary is ~ 5 x spin 
diffusion length. Also the grid spacing is varied until con- 
vergence is achieved. As an additional check we use the 
converged grid parameters for ID geometries to ensure 
that the theoretical results arc recovered. 

We first consider a three-terminal, NiFe (FM) and Cu 
(N), device with the following transport parameters 9 : 
cr NlPc = 8.62 x 10 6 /fhn, cr Cu = 59.52 x 10 6 /^m, A NlFc = 10 
nm, A Cu = 140 nm, a NiFo = 0.4, a Cu = and r c = flm 2 . 
For such a device the spin current is injected using the 
ferromagnet and detected in the normal metal. For this 
device we measure the amount of spin current injected 
into the normal metal as a function of the bridging dis- 
tance (d) between the ferromagnet and the normal metal. 
According to the ID SDD model the spin injection ratio 
(7) is given by, 



a 



7 



Ajy(T F (l-a 2 ) . ^ : 



(14) 



which does not depend on the bridging distance. Plot in 
Fig. [U shows 7 as a function of the bridging distance for 
a more general 3D device (Fig. [JJ inset). It can be clearly 
seen that the ID model is recovered in the limiting case 
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FIG. 2. (Color Online) Schematic representation of the various device geometries. Fig (a) and (b) represent the device geometry 
with varying contact area and thickness (H). Fig (c) represents the device geometry used to study the pinhole defects in 
tunnelling barriers. 



5£ 

^ 3 L 



•-•H= 


= 100 nm 




=500 nm 


m-m h= 


= 1000 nm 




=3000 nm 




=5000 nm 




0.5 



1 1.5 2 

5 2 

Contact Area (10 nm ) 



2.5 




Height ofSC layer (10 nm) 



FIG. 3. (Color Online) Spin injection ratio as a function of 
the contact area for direct contact device. Different curves 
are for different semiconductor heights. 



FIG. 4. (Color Online) Spin injection ratio as a function of 
semiconductor height for a direct contact device. Different 
curves are for different contact areas. 



where d > Xn (square, red-online). As d decreases 7 
increases rapidly and thus in order to achieve a better 
spin injection ratio the ferromagnet should be connected 
close to the normal metal. It should also be noted that 
in the ID model the only device dimension is along the 
flow of current and hence the ID model is quite limited. 



III. RESULTS AND DISCUSSIONS 

In order to show the potential of 3D modelling we 
study the effect of height (H) of the semiconductor 
and contact area (CA) [see Fig. [5] (a)] on n-GaAs with 
the following transport parameters: cr n _ G aA B = 10 5 /£lm, 
A n _GaA S = 1 Mm, Q^n-GaAs = and r c — fhn 2 . Kumar 
et ali£ have shown the effect of nano-pillar ferromagnet 
on the spin injection ratio by incorporating the effects 
of spreading resistance in the ID model, they could not 
show the effect of SC height since the model was es- 
sentially ID. It should also be noted that they used a 
contact area of ~ 12.5 nm 2 which is extremely difficult 
to reproduce experimentally and hence we study the ef- 



fect of varying contact area on the spin injection ratio 
(Fig. [3]). It can be seen from Fig. [3] that as the con- 
tact area increases 7 decreases rapidly. Due to the rapid 
decay a contact area of at least ~ 10 3 nm 2 would be 
required in order to achieve a direct contact spin signal 
into the device. According to our knowledge there has 
been no experimental evidence of direct contact spin in- 
jection into a semiconductor because most experimental 
geometries have a contact area of ~ 10 5 nm 2 . We also 
study the effect of sample height on the spin injection 
ratio (Fig. [4]). For a given contact area the spin injection 
ratio increases nearly linearly for small sample thickness, 
but as the sample thickness approaches the order of A sc 
the rate of change of the spin injection ratio reduces dras- 
tically. Thus in order to achieve better spin injection in 
direct contact devices a smaller contact area and a sam- 
ple height 3> A sc are required. Since extremely small 
contact areas are not feasible experimentally we inves- 
tigate the role of tunnelling barrier (AlO) and a thin 
metal (Cu) layer insertion for device geometry shown in 
Fig. [3 (b) . The parameters for the AlO tunnelling bar- 
rier are>22 r c = 1CT 7 Jim 2 and AS/E = 0.3 and the 
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FIG. 5. (Color Online) Spin injection ratio as a function of 
the contact area for AlO tunnelling barriers and Cu buffer 
layer for various semiconductor heights. 



FIG. 6. (Color Online) Spin injection ratio as a function of the 
pinhole coverage area. Square (red online) shows the result 
for a perfect tunnelling barrier without pinholes. The total 
contact area for this device is 0.5/im x0.5^m. 



thickness of the Cu insertion layer used is 50 nm. Tun- 
nelling barriers are an excellent solution to the conduc- 
tivity mismatch problem, but one of the main practical 
problems with tunnelling barriers is pinhole defect o 23 ' 24 . 
If we consider "ideal" tunnelling barriers without any 
pinholes then Fig. [5] shows that they are robust to vari- 
ations in contact area and N/SC height. But for small 
contact areas we can see that thin Cu films (~ 50 nm) 25 
can be excellent substitutes for tunnelling barriers. Here 
7 is measured at the Cu-SC interface, taking into ac- 
count the spin relaxation within the Cu buffer layer. It 
should be noted here that although using thin Cu film 
will decrease the effective spin diffusion length of the de- 
vice, thin metal films don't have defect problems like tun- 
nelling barriers. Hence only for small contact areas and 
where the spin diffusion length is not very important for 
the device thin metal films can act as excellent injectors 
of spin into the device. 

Lastly we discuss pinholes in tunnelling barriers and 
its effect on the spin injection ratio. Although some dis- 
cussions about pinholes^ have been made in the tun- 
nelling magneto-resistance (TMR) experiments based on 
a simple resistor model by Oliver et alj2£, there has 
been no work on this topic from the view point of the 
SDD model which is in general applicable to all spin- 
valve devices. In order to discuss pinholes we study a 
ferromagnet-semiconductor (FeNi-GaAs) interface with 
a single pinhole in the tunnelling barrier (AlO) as shown 
in Fig. [5] (c). The transport parameters used for GaAs 
art£: cr GaAa = 10 3 /f2m, A GaAe = 1 /mi, a GaAa = 0. Fig. [5] 
shows the effect of pinhole coverage area, i.e., the ratio 
of the pinhole area to the total contact area, on the spin 
injection ratio for the device. There are two competing 
effects here that determine the spin injection ratio, one 
of the conduction electrons passing through the pinhole 
(reducing 7) and second due to the tunnelling electrons 
passing via the tunnelling barrier (increasing 7). In or- 
der to understand the results obtained via simulations 



and thus the effect of pinholes, let us consider two chan- 
nels one for the conduction electrons (pinhole channel) 
and the second for the tunnelling electrons (tunnelling 
barrier channel). By assuming that these channels are 
independent it can be shown that the spin injection ra- 
tio of the device comprising of these two channels can be 
given by, 

jl.t Js.t _ Jl ,1 _ J2.4- Jl J2 

7 = : : : — = 7i^f Y~ + I2 -z — • 

(15) 

where ji, 2 , T ,4., Jx,a = ji.a.t + ji.a,^ are the up/down spin 
currents and total currents for the two channels and 
7i.2 = (ji.a.t - ji,2,J/(ji,2.T + ji.2,4.) are the spin injec- 
tion ratios of the two channels. 

Now for the case of the smallest pinhole if we consider 
the effect of area as shown previously (Fig. |3]) we obtain 
the 71 = 7.38%. Using the result for perfect tunnelling 
barrier (72 = 11.85%) and the fact that the pinhole acts 
like a short circuit causing most of the current to pass 
through the pinhole (Ji/(Ji + J 2 ) = 95% from simula- 
tions) we get 7 = 7.6% (from Eq. (fTof ). which is quite 
close to the value 7.26% obtained via simulations. The 
discrepancy of 0.34% seen is due to the fact that the tun- 
nelling barrier region now has a pinhole defect in it which 
was not considered while calculating j 2 . This result can- 
not be understood using the ID SDD model, which gives 
7 = 0.6%, since the effect of area cannot be taken into ac- 
count. Thus overall effect of pinholes is not just the sum 
of two individual ID channels. We also simulate more 
than one pinhole to see if there is any correlation between 
the different pinholes. We observe no such correlations 
at distances of ~ 300 nm, ~ 500 nm and ~ 800 nm which 
are typically the experimentally observable distances for 
such spin valve devices. Hence we conclude that the spin 
injection ratio depends only on the effective coverage area 
of the pinholes and not the number of pinholes present 
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in the tunnelling barrier. 



IV. CONCLUSIONS 

In summary we present the results of a 3D SDD model 
and show that the ID model fails to describe even the 
simple three terminal devices. We discuss the effects of 
N/SC height and contact area on the spin injection ratio 
of n-GaAs and show that direct contact spin injection 
is possible only for extremely small contact areas and 
height 3> A sc . We discuss the role of tunnelling barri- 
ers and show that thin metal films could be efficient spin 
injectors for small contact areas. Lastly the role of pin- 
holes is discussed and we show that the spin injection 
ratio depends only on the effective area of the pinholes 
and no correlation between the number of pinholes and 
7 is observed. 
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